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ABSTRACT 

Context. M87 is the first extragalactic source detected in the TeV range that is not a blazar. With the increasing performances 
of ground-based Cerenkov telescopes, we can now probe the variability in the y-ray flux at small timescales, thus putting strong 
constraints on the size of the emitting zone. The large scale jet of M 87 is misaligned with respect to the line of sight. A modification 
of standard emission models of TeV biazars appears necessary to account for the y-ray observations despite this misalignment. 
Aims. We explain TeV y-ray spectra and fast variability of M 87 by invoking an emission zone close to the central supermassive black 
hole, which is filled with several plasma blobs moving in the large opening angle of the jet formation zone. 

Methods. We develop a new multi-blob synchrotron self-Compton (SSC) model with emitting blobs set on a cap beyond the Alfven 
surface in the jet, at a distance of ~ 100 from the central engine to interpret the high energies inferred by new TeV observations. We 
present a SSC model that is explicitly adapted to deal with large viewing angles and moderate values of the Lorentz factor inferred 
from (general relativistic) magnetohydrodynamic models of jet formation. 

Results. This scenario can account for the recent y-ray observations of M 87 made by the High Energy Stereoscopic System (H.E.S.S.) 
telescope array. We find individual blob radii of about lO'"* cm, which is compatible with the variability on timescales of days recently 
reported by the H.E.S.S. collaboration and is of the order of the black hole gravitational radius. Predictions of the very high energy 
emission for three other sources with extended optical or X-ray jet which could be misaligned biazars still with moderate beaming are 
presented for one Seyfert 2 radiogalaxy, namely Cen A, one peculiar BLLac, PKS 0521-36, and one quasar, 3C 273. 

Key words, galaxies: active - galaxies: individual (M 87, 3C 273, Cen A, PKS 0521-36) - gamma rays: theory - radiation mecha- 
nisms: non-thermal 



C ■ 1. Introduction Due to the presence of a SMBH in the core and the pres- 

ence of the je t, M87 was deemed a n interesting candidate for 

- . M87 is a well-known nearby giant elliptical galaxy (z=0.00436, ^eV emission. EJBdiec et all (EUol reported an upper Umit ob- 

>^ ; I Smith etal. || 2000|) close to the center of the Virgo cluster, which ^^^^^ ^j^ij Whipple in 2000 and 2001, simultaneously with X- 

^ shows a multi-spectral jet, signature of an active galactic nucleus ^.^y fl^^g, observed by RXTE. HEGRA observed M 87 in 1998 

^ (AGN). Its jet IS one of the best known, at all scales, thanks to 1999 for a total exposure of 77 h aft er data quality selec- 

. P_ , Its nearby location and its strong synchrotron radiation in the ^^^^ dAharonian et alJlM:lBeilicke et alJl20 04). AA.lcr signif- 

optical band. M 87 is classified as FRI based on its radio mor- j^^^^g recorded and an integrated flux (E > 730 GeV) of 

phology. iWilson & Yana (120021) observed the jet with Chandra 33% Qiah was measured 

on July 29 and 30, 2000 and detected it up to a distance of ' Recendv. ^to^ian et all (l2006h have observed M 87 with 
-21 from the core in the X-ray band, which imphes that the jj^^ g^^^gy Stereoscopic System (H.E.S.S.flbetween 2003 
let IS not as strongly aligned along the Ime of sight (see also 2OO6 in 89 h live-time with a I3(r detection and discovered 
[ Reynolds et al. || 199 61) as in the case of biazars. variations on timescales of about 2 days, 10 times faster than 
At radio wavelengths, an impressive jet, which extends up t^at observed in any other waveband. This shows that the emis- 
to a few tens of kiloparsecs, can be seen. The central engine is gion j-ggion is very compact, with a dimension of the order of 
thought to be a supermassive black hole (SMBH) with a mass of ^ few Schwarzschild radii. Thes e observa tions, thus confirming 
Mbh ~ 3 X IO^Mq ( |Macchetto etal. || 1997|) . The scale length is t^e detection by HEGRA (Beilic ke et al.ll2005l) . are particularly 
thus Kg = GMbh/c ~ 4.5 X 10 cm - 1.4 x 10 pc. Using the interesting since M 87 is the first non-BL Lac extragalactic ob- 
Hubble Space Telescope (HST), [Biretta et al.| (1999) observed jg^t ever observed at TeV energy. Radio-loud galaxies contain 
superluminal apparent motions of about 4c-6c beyond 400 pc ^GNs with jets like biazars, but the jet emission is less strongly 
for the internal knots, between 1994 and 1998, thus confirming boosted due to larger viewing angles between the jet and the 
that the jet is relativistic. They conclude that the jet is oriented observer's line of sight. It is therefore a challenge for standard 
within 19° of the line of sight. models of TeV biazars to explain the very high energy (VHE) 
emission of M 87. 
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In this paper, we present a modified synchrotron self- 
Compton (SSC) scenario to explain the VHE emission of 
M87. Classic SSC models (e.g., G ould 1979; Inoue & Takahara 
1996; Bloom & Marsche3 ll996riChiaberge & Ghisemnilll999F 
Katarzvhski et al.. ,200 lb are applied to blazars, which are 



beamed sources, and cannot account for the observations of ra- 
diogalaxies like M87. Our goal is to further develop one of 
these models to reconcile beamed and unbeamed sources in the 
same framework of models. Such propositions for unification of 
AG Ns have already been studied cons idering orientation effects 
(e.g. lAntonuccil 119931; lUrrv & Padova ni 1995), or radio/X-ray 
power among BLLac obje cts, flat-spectrum radio-loud quasars 
(FSRQs) and FRIs (e.g . iFossati et al.1 119981: iGhisenini et all 
119981: ICapettiet al.ll2000h . 

A short description of the leptonic blob-in-jet model for TeV 
blazars is found in Sect. [21 and its development and application 
to M 87 are described in Sect. [3] and Sect. |4] In the framework 
of misaligned BLLac-like objects, we then try to predict VHE 
fluxes for objects with optical/X-ray extended jets in Sect. |5] 
Implications on unification schemes of AGNs are discussed in 
Sect.|6l 



2. "Blob-in-jet" leptonic SSC model 

We intend to model the multiwavelength spectrum of M 87 in 
the framework of a quasi-homogeneous SSC scenario, success- 
fully used to account for the overall emission of blazars, such 
as Mrk 501 and Mrk 421. Our model relies on the basic sce- 
nario presented in Katarzvnski et al. (2001, 2003, and references 
therein) who give the details of the computation of the radia- 
tive transfer and emission by SSC processes in a single spher- 
ical blob of plasma moving at relativistic speed along the jet 
axis. The blob, immersed in a uniform magnetic field, is as- 
sumed to be located inside the jet, close to the central engine. 
An inhomogeneous conical extended jet model is also used to 
explain the emission from radio to u ltraviolet wavelengths (see 
Sect. 2.2 in iKatarzvhski et alJl200lL for more details). The ab- 
sorption by the infrared extragalactic background light at VHE 
is taken i nto account and mod eled using the estimations as de- 
scribed in lStecker et alj (l2006l) and references therein. Here we 
model only nearby active galactic nuclei (AGNs) and hence this 
effect can be neglected. The blob-in-jet model is particularly 
well adapted to the description of blazars, for which the jet is 
very close to the line of sight. In the following, the assumed 
cosmology is Hq - 70 kms^^Mpc ' for an Einstein-de Sitter 
universe, with Qa - 0.7 and Q,„ - 0.3. 

We assume that the population of electrons, which is re- 
sponsible for the non-thermal emission in leptonic models, has a 
number density that can be described by a broken power-law: 



Kiy-"' 



< r < Tc 



[cm-3] 



(1) 



where K2 - Kij'^^^"' and 7 = E/mc^, where m is the electron 
mass and E its energy. These electrons radiate up to the X-ray 
range through the synchrotron process, and then re-interact with 
their own emitted photons by inverse Compton (IC) scattering, 
which is the so-called synchrotron self-Compton process. This 
synchrotron emission comes from a population of electrons dif- 
ferent from those producing the radio-IR emission of the ex- 
tended jet. 

The SSC model has 8 significant parameters. The macro- 
physics processes are described by the magnetic field B, the 
radius of the emitting blob r/, and the Doppler factor St - 



[Th{l -ySfc cos 6*)]"', where /^i, is the speed of the moving blob 
in c unit, F;, = (1 - y6^)"'^^ is the blob Lorentz factor and is 
the viewing angle. The radiative processes are parametrized by 
the description of the population of emitting particles, with the 
parameters K[, ybr, 7c, "i and n2 from Eq. ([TJ. The value of ymi,, 
is not crucial for the interpretation of the spectral energy distri- 
bution (SED), nor is jc, although it can become very relevant in 
cases where the X-rays have a hard slope with a spectral diff'er- 
ential index o- < 1 (in the common fy oc v^" notation). All these 
parameters can be constrained when detailed spectral data are 
available for a wide frequency range. 

In the present case, the spectral coverage of the nucleus of 
M 87 is sparse and we need to find other ways to constrain the 
parameters. One important constraint comes from the variability: 



rb < 



c6b 
1+z' 



Af, 



obs 



(2) 



where Afobs is the variability timescale in the observer frame, 
implying rb/di, ^ 5 x 10'^ cm for M 87. 

The region of emission is then assumed to be close to the 
SMBH, to fulfill the variability constraint within magnetohy - 
drodynamic (MHD) jet models. For instance. iMcKinnevI (l2006h 
models the jet formation zone using general relativistic mag- 
netohydrodynamic simulations, applicable to GRBs, AGNs as 
M 87 and black hole X-ray binaries. He describes the broaden- 
ing zone of the jet in the vicinity of the central black hole, and 
finds the Alfven surface at ~5Q rg. We assume that the emission 
zone is located slightly above this surface to allow shocks and 
Fermi acce leration processes to develop in the jet. The results of 
iMcKinnevI (|2006) further constrain some of our parameters for 
a distance of ~100ri, from the SMBH, such as the value of the 
Lorentz factor F^ $ 10 of the plasma blobs, the magnetic field B, 
and the half-opening angle (p{r) of the jet given by his Eq. (24). 

The case of M 87 is of particular in terest since its jet is ex- 
ceptionally well mapped in radio VLBI. lBiretta et all (|2002|) ob- 
served the core of M87 in VLBI in February 1995 and March 
1999, and showed that the opening angle increases quickly with 
decreasing distance to the core region, at the 0.01 pc (~lQrg) 
scale. Such a widening at the base of the j et was also observed in 
CenA at the 0.1 pc (-19000 r^) scale by iHoriuchi et al.l(l2006h 
using the VLBI Space Observatory Programme. A broaden- 
ing zone in the jet formation region is found as well in MHD 
simulations. Moreover, the recently detected short term TeV 
variability seems to exclude the extended and outer regions as 
the source for VHE emissio n in M87. This is also argued by 
Georganopoulos et al.' ("2005) who present a modified leptonic 
model applied to M 87 which takes into account a deceleration 
of the inner flow along the base of the jet. So it appears quite 
natural to assume that the VHE y-rays are emitted in the core 
widened jet region. 

We can then imagine that there are blobs of plasma, harbor- 
ing very high energy electrons and propagating in the widened 
jet formation zone, that are dragged along with the bulk jet out- 
flow. In the case of misaligned objects such as M 87, this can 
easily result in one blob moving along the line of sight and thus 
having about the same Lorentz factor as for blazars, allowing 
to reproduce the TeV emission in the framework of classic SSC 
models. However a model with a single relativistic blob moving 
and emitting exactly towards the observer would be statistically 
unlikely. 
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Fig. 1. Geometric side view of the jet formation zone ("on-blob" 
case). Vj represents the velocity of the jet, the blobs are num- 
bered from to 6, vo^e are the velocities of the blobs, 6 is the 
viewing angle with respect to the jet axis, (p(r) is the opening an- 
gle, rh is the radius of an individual blob and R^ap is the distance 
of the blobs from the SMBH. In this sketch, the contributions to 
the total flux of the blobs n°2 and 3 would be rejected, since they 
lie outside the jet. 



3. Multi-blob model 

A way to deal with this statistical issue is to assume that the 
emission zone is a spherical cap centered on the SMBH, limited 
by the sheath of the jet and filled with sever al similar homo- 
geneo us blobs. Consequently, as mentioned in lAharonian et al] 
(120061) . we can consider differential Doppler boosting in the jet 
formation zone, near the core region. 

This cap is located at a given distance T^cap from the SMBH, 
which is a new free parameter in our mo del. However, R^ ^^ can 
be constrained by MHD simulations (e.g. McKinneyl2006l) if we 
assume that it is located slightly above the Alfven surface, which 
is at about 50-100 from the SMBH. This surface is continu- 
ous, but does not need to be homogeneous. We model it with 
a pattern of several blobs, whose individual radii are typically 
smaller than in the case of the "blob-in-jet" scenario. 

For the sake of simplicity, we assume in this zone the pres- 
ence of 7 blobs, one central blob and 6 further blobs distributed 
on a hexagon, located at 1 00 from the SMBH, with macro- 
scopic parameters derived from iMcKinnevI ( l2006l) as specified 
in Sect.|2] This choice for the number of blobs is justified by the 
fact that the resulting diameter of the cap is of the same order 
of magnitude as the characteristic size of the emitting zone from 
previous studies. 

The smaller the radius of the individual blob, the more our 
model resembles a continuous zone model. The choice of dis- 
crete adjacent blobs leads to two extreme situations: 

- "Inter-blob" case: the line of sight passes exactly through the 
gap between three blobs. 

- "On-blob" case: the line of sight is exactly aligned with the 
velocity vector of the central blob, called n°0. 

The Doppler factor for each blob n°i (where / e |[0; 6J) is 
now defined as: 



3 



Fbil- /3b cos ad 

where a, is the angle between the velocity vector of blob n°/ 
and the line of sight. If the line of sight is between three blobs 
("inter-blob" case), then these blobs have the same Doppler fac- 
tor and their contribution to the total flux is equal, while the 
4 other blobs have contributions to the total flux that decrease 
with increasing blob radius rt- If the line of sight is aligned 
with the velocity vector of the central blob ("on-blob" case), 
then the highest Doppler factor is 5". In that case the six other 
blobs have all the same Doppler factor, and although smaller 
than 5" their individual contributions are not negligible in the 
total ol3served flux, especially if the seven blobs are all moving 
along in the same directior0 It should be noted that some mod- 
els involve acceleration as the je t is collimated (e.g. Melia e t aTl 
120021 IVlahakis & Konigllllool . that is at the parsec scale. In 
such models, a gradient in flow velocity across the width of the 
jet can also be present, but is usually small compared to the ra- 
dial velocity profile. Therefore we choose here to neglect this 
transverse gradient, which is of second order for our purpose, 
and we assume that all the blobs have the same Lorentz factor, 
although they are ejected at slightly diff'erent angles. 

Figure [1] describes the geometry of our model in the "on- 
blob" case. The central blob is moving along the line of sight, 
and the six other blobs are each moving along a direction slightly 
different from the blob n°0 by an angle da. The angle da is given 
by c/a = 2 arcsin (rfc/Z^cap)- The viewing angle is defined as the 
angle between the line of sight and the jet axis for this multi-blob 
model. The individual radius assumed equal for all blobs is a 
free parameter of our model. Depending on the observed angle 
6, it can happen that in the simulation a blob moves outside the 
jet and is therefore neglected. 

We can then compute the radiativ e transfer of each blob i n 
its own source frame, as explained in iKatarzvnski et aTl ( 1200 lb . 
For each blob, the seed photons for the inverse Compton scat- 
tering are those generated within the blob from synchrotron ra- 
diation. The total flux in the observer frame is the sum of the 
contribution of each blob. We neglect the contribution of the 
eventual blobs for which 6b < I, which is possible if a, > 
arccos[(r,, - l)/(ri/?,,)]. 

Figure |2] shows the SED of the contribution of the single 
central blob in the "on-blob" case (thin solid line), while the 
thin dashed line shows the contribution of the six adjacent blobs. 
Adjacent blobs have the same Doppler factor, so they all con- 
tribute equally to the SED. The bold solid line shows the total 
flux of the system, which is the sum of the contributions of the 
seven blobs. The bold dashed line shows the sum of all the con- 
tributions in the "inter-blob" case, where the line of sight is ex- 
actly between three blobs. If r/, is sufficiently large, the closest 
blob to the line of sight in the "on-blob" case completely dom- 
inates the apparent flux. This can be the case in almost all the 
situations we study here. 

4. Application to M87 

4.1. The observed SED 

To construct the SED of the core jet, we carefully selected the 
following data from the literatu re. We use y-ray obse rvations 
of 2004 and 2005 by H.E.S.S. ( lAharonian et al.ll2006l) . shown 
in black and gray points respectively in the next figures. The 

- which is the case for rj, < rg. 
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loR(r) (Hz) TesL 

Fig. 2. Illustrative example of the SED of the multi-blob model. 
The thin solid line shows the contribution of the central blob, 
the most strongly beamed blob. The thin dashed line shows the 
contributions of each of the six other blobs, and the bold solid 
line shows the sum of the contributions of all the blobs in the 
"on-blob" case. The bold dashed line shows the sum of all the 
contributions in the "inter-blob" case. The closest blob to the 
line of sight (thin solid line) is overwhelmingly dominant for 
sufficiently large values of rt,. 



Whipple upper li mit at 400 GeV observ ed between 2000 and 
2003 is taken from lLe Bohec et al.l(l2004 . The H EGRA point at 
730 Ge V obtained in 1998 and 1999 is taken from lBeihcke et all 
(12004 . 

The Chandra data from the nucle us region obtained on April 
20 and July 30, 2000 are taken from lPerlman & WilsonI (l2005h 
who conclude that the nuclear X-ray emission originates from 
the jet and could extract a 1 arcsecond nucleus spectrum thus ex- 
cluding HST-l . The XMM-Ne wton data taken on June 19, 2000 
were found in lBohringer et al.l (1200 lb . It should be noted that the 
XMM-Newton data have the same spectral slope, but a higher 
flux density than the Chandra d ata, and they do inclu de HST- 
1 within a 4" extraction region. iPerlman et al.l (1200 lb provide 
data taken by the Gemini North telescope at 10 microns in May 
2001. The observations by the HST in the o ptical and UV bands 
in 1991 are reported in ISparks et al.l (Il996i) . The observations 
of the core of M 87 by the VLA in the radio band and by the 
Palomar observat ory in the optica l band between 1979 and 1985 
were found in Biretta et al.l (|1991|) . 

We also take into ac count upper limits in y- ray by EGRET 
betwee n 1991 and 199 3 ("Sreek umar et al.ll 19941) . and in UV by 
EUVE dBerghofer et al . 2000). All the other data are taken from 
the NASA/IPAC Exti-agalactic Database (NED). 

Although we do not have simultaneous data, the Chandra 
data of 2000 and the H.E.S.S. y-ray data of 2004 both correspond 
to low states of activity in their own spectral range. We therefore 
associate them, assuming that they are representative of a typ- 
ical low state. Indeed, regarding the X-ray data, the Chandra 
observations of 2000 correspond also to the lower state of ac- 




108(1-) (Hz) M 87 

Fig. 3. Tentative modeling of the SED of M 87 within a stan- 
dard blob-in-jet scenario, with (5^ = 3.07 presented by the dashed 
green line. The solid blue line shows the SED of M 87 emitted 
by a single blob moving along the line of sight in the jet forma- 
tion region, with 6b - 8, which can describe the data. The data 
points that bring direct constraints to our model are shown in 
black, the other less constraining points are represented in gray. 
The black line in the optical band shows the host galaxy, as- 
sumed to be elliptic. This co ntribution was computed using re - 
sults from the code PEGASE dPioc & Rocca-Volmerangell997h . 
The radio data shown by gray empty circles, obtained from the 
NED and with fluxes ~10"'^ergcm"^ s"', come from the ex- 
tended kiloparsec-scale jet and radio lobes. The black line from 
radi o to UV/X represents a mode l of the extended inner jet 
(see iKatarzvnski et al.ll200lL l2003h . with the corresponding ra- 
dio data with fluxes ~10"^^ ergcm"^ s ' reported in black. The 
bump peaking at -lO'^^Hz is due to synchrotron emission and 
the bump peaking at ~10^^ Hz is due to inverse Compton, both 
from the VHE zone (see columns 1 and 2 of Table[T]for the corre- 
sponding parameters). [See the electronic edition of the Journal 
for a color version of this figure.] 

tivity recently published (iPerknan et al.ll2003l) . For the models, 
we choose to take into account the mean spectral slope of the 
Chandra data. The radio to optical/UV data are also not simul- 
taneous with the y-ray data, but this is not problematic since the 
radio contribution is thought to come from the extended jet, with 
characteristics different from the VHE emitting zone. 

4.2. The SSC models 

In the case of M 87, the observati on angle 6 between the jet axis 
and the line of sight is at most 19° dBiretta et all 19991) . The blob- 
in-jet model cannot describe correctly the VHE emission of the 
source, as it would require very high Doppler factor, which is not 
sustainable. This is illustrated in Fig. |3] The dashed green line 
presents the best solution for the SED of M 87 within the blob- 
in-jet scenario described in Sect. |2] with 61, = 3.07 assuming 6 = 
19° (see column 1 in Table[T]for the corresponding parameters). 
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Fig. 4. SED of M 87 within the multi-blob scenario for small 
blob radius (ri, = 1.5 x lO'^'cm, see column 3 in Table [U. The 
two extreme "on-blob" and "inter-blob" cases are identical here 
because of the small value of r;,. [See the electronic edition of 
the Journal for a color version of this figure.] 



Greater values of 6h are not allowed here because of the large 
value of 0. We can obviously see that this model cannot account 
for the VHE emission. The well defined X-ray slope deduced 
from the observations strongly constrains the second index n2 
in the electron energy distribution (see Eq. ([T|i) and significantly 
reduces the parameter space. 

The solid blue line in Fig.[3]presents the resulting SED emit- 
ted when one consider a single blob moving along the Une of 
sight in the jet formation zone with 6h - 8. The correspond- 
ing parameters can be found in column 2 of Table [U where is 
defined as the angle between the line of sight and the velocity 
vector of the single blob. Obviously this model describes the ob- 
servations much better However, as pointed out in Sect. |2l it is 
based on an ad hoc assumption. Moreover, it is difficult to "keep" 
the generated IC bump below the EGRET upper limit although 
we assume a low state for the activity of the AGN. 

Throughout this paper, our results are not fits to the data, but 
rather solutions of models which are meant to describe best the 
data. Our purpose is to figure out whether our model can describe 
correctly the current available data for different objects. We do 
not intend to fine-tune the parameters of our model but to sort 
orders of magnitude out for these parameters. 

One SED of M 87 generated within the multi-blob model is 
presented in Fig. |4] with parameters very similar to the single 
blob model of Fig. [3] (see column 3 in Table [TJ. Since the for- 
mer is a generalization of the latter, the resulting spectrum is 
rather similar, as one would expect. In this case, the value of the 
individual blob radius r/, is so small that all the blobs are mov- 
ing close to the line of sight. The "on-blob" and the "inter-blob" 
cases give the same contribution to the SED and are overlaid in 

Fig. a 

The blob radius is rather small in this case, resulting 
in a VHE emitting zone smaller than the Schwarzschild ra- 
dius. It should be noted that features small compared to the 
Schwarzschild radius, possibly responsible for the VHE emis- 
sion, can develop beyond the Alfven surface due to turbulence 
or reconfined shocks, but this issue is beyond the scope of this 
work. 
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Fig. 5. SED of M 87 within the multi-blob scenario, assuming a 
low magnetic field {B - 10 mG, see column 4 in Table[T]i in blue 
lines. The solid lines show the "on-blob" case while the dashed 
lines represent the "inter-blob" case. The green lines show a 
solution for the high state observed by H.E.S.S. in 2005, with 
the same parameters as for the observations of 2004, except for 
rb = 8.0 X 10'^ cm, Ki = 2.2 x lO^^cm^^ and n2 = 2.5. [See 
the electronic edition of the Journal for a color version of this 
figure.] 



However it is commonly believed that the size of the VHE 
emitting zone cannot be much smaller than the Schwarzschild 
radius which is a natural scale for the processes in the vicinity of 
the SMBH. Moreover the emitting zone must be large enough to 
allow the acceleration of particles to develop. Very small blobs 
may disappear rapidly, in ~ 10 minutes due to adiabatic expan- 
sion which is especially important in the broadened zone of the 
jet. However a long, stable emission is possible, even from small 
blobs. The emitting zone can be located at a stable stationary 
shock front, above the Alfven surface. It initiates the accelera- 
tion and thus the radiation of particles of a large number of small 
blobs continuously crossing the shock, thus providing a quies- 
cent background of VHE emission. Density fluctuations in the 
injection of material could then generate flares as seen at VHE. 
In fact, the only problem with small blobs is that in this case the 
paving of the jet is not complete because of the discretization 
applied in our code. 

In order to be more conservative and to fulfill the constraint 
rb ^ rg, we analyze another possibility with a low magnetic field. 
It is presented in Fig.|5]in blue lines with associated parameters 
in column 4 of Table[T] In that case, the result predicted by MHD 
models with a strong magnetic field in the vicinity of the cen- 
tral engine is not strictly fulfilled. However a local decrease of 
the magnetic field can be achieved by a simple expansion of the 
emitting zone. This solution may appear preferable for conser- 
vative reasons on the size of the emitting zone. 

A satisfactory solution for the high state observed by 
H.E.S.S. in 2005 is also possible in this case within the multi- 
blob scenario as shown in Fig. |5] (green lines). Interestingly we 
predict a radical change of the X-ray regime. Clearly the param- 
eters are not well constrained here due to the lack of any simulta- 
neous data, especially in the X-rays. However this illustrates the 
capability of the multi-blob model to generate spectra that are 
sufficiently hard in the TeV range to reproduce the most recent 
H.E.S.S. data. A slight modification of standard SSC models for 
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Table 1. Parameters used in the different generated SEDs. 
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TeV blazars appears therefore successful and can account for all 
VHE data on the radiogalaxy M 87 that are available up to now. 

5. Predictions for other radiogalaxies with 
optical/X-ray extended jets 

We now apply the SSC multi-blob model to different sources 
that have the peculiarity to show an extended optical or X-ray jet, 
which suggests at least a moderate beaming towards the observer 
as in the case of M 87. This allows to confront the multi-blob 
scenario to other types of AGNs and to predict whether these 
sources can be detectable at VHE or not by present Cerenkov ar- 
rays. We choose three AGNs not belonging to the genuine blazar 
class, presented here with increasing viewing angles, such that 
their fluxes are less and less boosted by relativistic effects. We 
stress that not all the data presented here are simultaneous and 
that the sources undergo large variations. 

5.1. 3C273 

3C 273 (z=0.158. [s"trauss et al.lll992h is the first qu asar that was 
identified as a high-redshift object ('Schmidtj [r963h and the best 
studied. It hosts a SMBH whose mass is at le ast ~2.0xlO'^Mp) 
as in ferred from studies of Balmer lines (iPaltani & Tiirleii 
|200 5|) and its maxim um acceptable viewing angle is about 
15° (llJnwin et al.lll985h . 

Blazars display featureless X-ray contribution but radio- 
galaxies ca n have a much mor e complex environment at low en- 
ergies (e.g. lGrandi et al.l[2006h . However since the purpose here 
is to model the non-thermal contributions of these objects, we 
decided to consider the X-ray contribution as dominated by the 
jet emission, keeping only a feature in soft X-r ay as a signature 
of the accretion disk (Grandi & Palumboll2007b . 

Since this source is highly variable, one needs to be careful 
to select simultaneous data. The X-r ay data presented her e are 
BeppoSAX observation^ taken from lGiommi et all (|2002|) . We 
report in red in Fig.|6]the upper limit at 3cr obtained by H.E.S.S. 
in 2005 tAharonian et aTluOOSh. All the other dat a points pre- 
sented in gray are taken from iTiirler et al] d 19991) who report 
an average spectrum compiled from 30 years of observations. 
Observations by BeppoSAX on January 13, 15, 17 and 22, 1997 
lie within the observation period of CGRO taken by COMPTEL 
and EGRET between December 10, 1996 and January 28, 1997 

^ see |http : //www . asdc . asi . it/blazars/ 1 
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Fig. 6. SED of 3C 273 wifli anticipated VHE flux. The solid blue 
line shows the average of the "on-blob" and "inter-blob" cases 
(see column 5 in Table[T]for the corresponding parameters). The 
dashed line represents a simple blackbody model to illustrate the 
contribution of the big blue bump component. The upper limit 
obtained in 2005 with the H.E.S.S. telescope array is shown in 
red. The blue lower limit shows the expected CTA sensitivity in 
50 h of observation. [See the electronic edition of the Journal for 
a color version of this figure.] 



( Collmar et al.l l2000h which are also reported in iTiirler et"ari 
Jl999h . We have thus simultaneous data for the X/y-ray bump 
in one of the highest state, which puts an important constraint on 
the models. 

The nature of the X-ray emission of 3C 273 is still an unre- 
solved issue. The high frequency bump, thought to be probably 
due to IC emission, would present a peak at a rather low fre- 
quency compared to other AGNs. This implies that either 3C 273 
does not emit at detectable VHE levels, or that the nature of this 
bump is synchrotron, implying then the presence of a hypotheti- 
cal IC bump at ultra high energies. 

Assuming that the hard X-ray emission is indeed due to in- 
verse Compton emission, the solid blue line in Fig. |6] shows 
the SED of 3C273 with the predicted VHE flux (see column 
5 of Table [T] for the corresponding parameters). Modeling the 
SED of 3C 273, we note that changing the value of y^^^ dramati- 
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cally affects the X-ray flux by increasing it with decreasing ymin, 
and hence could explain some X-ray flares, as also suggested 
by Georganopoulos et al. ( 2006) in the case of external inverse 
Compton emission on the Cosmic Microwave Background. The 
value of 7min is strongly constrained here by the precise shape 
of the X-ray spectrum. In the following figures, the V-shaped 
curve at VHE (~10^^Hz) shows the H.E.S.S. sensitivity limit 
for a detection of 5cr in 50 h of observation for a source at a 
mean zenith angle of 30°. The expected sensitivity of the next 
generation CTA project of ~0. 1 % Crab flux at 1 Te V for 50 h of 
observation is shown as a blue lower limit. 

Since the results for the "on-blob" and "inter-blob" scenarii 
are not very different in this case, we show with a blue line the 
average of the two. We should stress that low frequency data 
are not simultaneous with the X/y-ray data that we selected. 
The modeled synchrotron bump has a higher flux density than 
the optical data because we are dealing with a high s tate of ac- 
tivity in y-rays as reported by IColknar et al.l ilQQu ). We only 
predict a very marginal detection of 3C273 by H.E.S.S. in its 
low energy range, depe nding on the energy threshold (but see 
also iGeorganopoulos et al. 2006). Furthermore it should be re- 
called that ICollmar etalT EoOO) report an active state in y-rays 
an d lTiirleretalJ (120061) a high level of the non-thermal emission 
at the epoch of the data we are considering. So even in a high 
state we do not expect a strong level of VHE y-ray within our 
scenario. 

A strong detection of 3C 273 at VHE with the current gener- 
ation of Cerenkov arrays would be difficult to explain within our 
SSC nuclear scenario. A possibility would be to invoke dispar- 
ity among the various emitting plasma blobs. Our model shows 
the presence of a well peaked inverse Compton bump; differ- 
ent magnetic fields or electron energy distributions among the 
blobs could result in a tail of the IC bump at VHE that could ac- 
count for a VHE detection. An alternative would be an extended 
emission due to external inverse Compton radiation, which is 
then expected to be not very variable. In all cases, further obser- 
vations with GLA5r( 10 keV-300GeV) and H.E.S.S. II, which 
will extend the spectral domain of H.E.S.S. I down to ~20GeV 
with a better sensitivity, are required to disentangle the different 
plausible scenarii. 

5.2. CenA 



Cen A is the nearest radiogalaxy (z=0.0018. lGrahamlll978h and 
one of the best studied. The presence of an AGN in CenA is 
evident from the radio-ba nd to the y-rays. Observations from 
CGRO ( Kinzer et al.lll995l) show a bump that seems to peak at 
~200keV as pointed out bv! Steinle et al.. (1998). We thus have a 
real constraint for the parameters of our model, particularly with 
regard to the description of the population of electrons. As in 
the case of 3C273, the I C bum p peaks at a rather low energy, 
leading Chiaberge et al.l (1200 lb to note that Cen A could be a 
misaligned low-energy peaked BL Lac (LBL) object. The value 
of the vi ewing angle of the j et is still a controversial issue. For 
instance, iTingav et al.' ('1998) claim d ~ 50°-80° for the parsec- 
scale jet, whereas Hardcastle et al. (2003) find 6 ~ 15° for the 
lOOpc scale jet . We choose here to take an intermediate value 
of 6 ~ 25° (See lHoriuchi et all (l2006h for a discussion about the 
different values for the viewing angle of Cen A found in the liter- 
ature). The SMBH mass inferred from gas kinematical analysis 
usin g a [S III] line is ~ 1 . 1 xlO^Mp) (Marconi et al. 2006), but see 
also lHaring-Neumaver et al.l(l2006l) who give Mbh ~ 6 x WMq 
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Fig. 7. SED of Cen A within the multi-blob scenario. Solid lines 
show the "on-blob" emission; the "inter-blob" cases are rep- 
resented in dashed lines. In blue we show a model assuming 
that the soft y-ray data are inverse Compton emission, while in 
green we assume a synchrotron emission to describe them (see 
columns 6 and 7 respectively in Table [T] for the corresponding 
parameters). [See the electronic edition of the Journal for a color 
version of this figure.] 



The data sample c hosen here is almost the same as in 
I Chiaberge et al.l (l200lh . The data were retrieved carefully to 
take only the nucleus into account. The strongl y constrain- 
ing CG/?(9/COMPTEL y-ray data are taken from Steinle et~an 
( 1998); Rothschild et al. (2006) provided RXTE and INTEGRAL 
data; //5r/NICM0S and WFPC2 data, which we have carefully 
dereddened, ai-e from Marconi et al. (20000; SCUBA data at 
800 jum are taken from Hawarden et al. ( 19 93|); ISO and SCUBA 
(450//m and 8 50 ^m) data are fro m Mirabel et alj (11999 ); VLA 
data are from lBurns et all (ll983l) . lEvans et all ( |2004 report X- 
ray observations by XMM-Newton on February 2 and 6, 2001 
and by Chandra on May 9 and 21, 2001 with a photon index 
r = 2 for the parsec-scale jet component. Data from the NED 
are shown as non-constraining points (in gray) for comparison. 
The H.E.S.S. upper limit based on observation s in 2004 with a 
live time exposure of 4.2 h is reported in red dAharonian et al.l 
2005). 

We should also point out that as Cen A harbors a strongly ab- 
sorbing dust lane, and because this source is extremely close and 
well resolved, the X-ray data should be then taken only as upper 
limits. The contribution of the nuclear jet might be contaminated 
by the accretion disk of the AGN and by the X-ray binaries that 
are resolved in this object. In this case we would have only poor 
constraints on the emission process. We assume here that all our 
selected high energy data come from SSC processes. 

Figure |7] shows the SED of Cen A applying the multi-blob 
model in two cases, (/) assuming that the y-ray peak observed 
by CGRO is inverse Compton radiation (solid and dashed blue 
lines) or (//) assuming it to be synchrotron (solid and short 
dashed-long dashed green lines). 

We should als o point out that the previous study by 
ISteinle et alj (Il998h reports a variability in soft y-ray of about 



using a [Fe II] line. 



• * The reddening co rrection was applied using the factors carried out 
bv lMarconi et al J ( 120001 themselves. 
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10 days, implying rt, < ly. lO'^cm (see Eq. which is well 
satisfied by our parameters. 

Given the results of our model in the first scenario with an IC 
bump in soft y-rays {blue lines), the SSC emission of the central 
region would definitely not provide a flux sufficiently high to be 
detectable at VHE (see column 6 in Table [T]i, at least for a SSC 
emission dominated by the nucleus. This holds even in the case 
of huge v ariations of the nuclear emission. This conclusion con- 
curs with lStawarz et al.l (l2003h . who do not expect SSC emission 
by the nucleus or by the base of the jet of Cen A, but do expect 
VHE emission that could be detectable by current imaging at- 
mospheric Cerenkov telescopes facilities in the case of an exter- 
nal inverse C ompton emission pro cess on the host galaxy pho- 
ton field. For lStawarz et al.l (l2006h . many Fanaroff'-Riley type I 
(PR 1) radiogalaxies like Cen A could be TeV sources for which 
the weak nuclear y-ray emission would be absorbed and re- 
processed by inverse Compton on the starlight photon field, thus 
generating an isotropic y-ray halo. In our model, the lack of si- 
multaneous data prevents us from further constraining the syn- 
chrotron bump, which has here a higher flux density than the 
selected data since we are considering a state of high y-ray ac- 
tivity. 

In solid ("on-blob" case) and short dashed-long dashed 
("inter-blob" case) green lines in Fig. |7] we present a SED of 
Cen A assuming now that the soft y-ray peak is of synchrotron 
origin. In this case (see column 7 in Table[T]i, we expect a detec- 
tion of the core of Cen A at VHE by the H.E.S.S. tele scope array 
within 50 h of observation. It should be noted that iBai & Leg 
(1200 ih also predicted the synchrotron bump to be in the soft y- 
ray range and the inverse Compton bump to peak around 1 TeV 
in the context of SSC models, which comforts our latter model. 

5.3. PKS 0521-36 

PKS 0521-3 6 is a FSRQ object with an optical jet located at 
z = 0.055 d liillilll The central SMBH h as a mass of 
~3.3x lO^Mr:^ (IWoo et al.l[2005h. lPian et aPdiggeh. and more re- 
cently iTingav & Edwards (12002"), mention the absence of super- 
luminal motions in its jet, contrary to the case of 3C273, im- 
plying that the beaming effect is much less important and thus 
strengthening their result on the viewing angle . Indeed the only 
constraint on the jet orientation comes from Plan et all d 19961) 
who deduce 6 ^ 30° + 6° from SSC models. We should also note 
that PKS 0521-36 seems to osci llate between a Seyfert-like and 
a BLLac state (e.g. lUlrichlll981h . making this source difficult to 
interpret within a pure non-thermal scenario, especially since we 
are confronted with non-simultaneous data. 

BeppoSAX obse rved PKS 0521 -36 on October 10, 1998 
(iGionuni et al.ll2002h (black points in FigM, and the SwiftPiRT 
measurements (green points) taken on May 26, 2005 were ob- 
tained through the Online Analysis TooQ The data points in 
gray are from the NED. We report in red the upper limit at 2cr ob- 
tained in 89 h from observatio ns by CANGAROO between 1993 
and 1996 dRoberts et aljl998i) . We fu rther used the EGRET data 
between 30MeV and 500 MeV from lHartman et alJ (Il999h and 
taken between July 12, 1994 and August 01, 1994. The blazar 
PKS 0521-36 is associated with ffie source 2EG J05 24-3630 in 
the Second EGRET Catalog dThompson et al.ll 19951) . but during 
cycle 4 this source was found to lie outsi de the 99% co nfidence 
contour of EGRET. However, like Fosch ini et alj d2006l) . we as- 
sume in this work the identificat ion with PKS 0521-36 to be 
vahd, which is also pointed out bv lTornikoski et alj d2002h . 
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PKS 0521-36 



see |http : //www . asdc . asi .Tt] 



Fig. 8. SED of PKS 052 1-36 with anticipated VHE flux mod- 
eled in the multi-blob framework. The solid blue line is the "on- 
blob" case while the dashed green line shows the "inter-blob" 
emission (see column 8 in Table[T]for the corresponding param- 
eters). The CANGAROO upper limit is reported in red. [See the 
electronic edition of the Journal for a color version of this fig- 
ure.] 

Figure [8] presents the SED of PKS 0521-36 with the antic- 
ipated VHE emission (see column 8 in Table [T] for the corre- 
sponding parameters). It seems unlikely that the X-rays are due 
to inverse Compton radiation since the inverse Compton scatter- 
ing would be with photons from the radio/optical contribution, 
which is not very variable, coming from an extended part of the 
jet. Since the X-rays show a high degree of variability, they most 
certainly come from a compact region and are of synchrotron 
origin. If the X-ray emission is truly due to the synchrotron pro- 
cess, we predict that this BLLac object should be marginally 
detectable by the present H.E.S.S. telescope array, and easily de- 
tectable by H.E.S.S. II and by the next generation of Cerenkov 
arrays, such as the CTA project, which will detect sources down 
to -0.1% of the Crab flux. Furthermore, if PKS 0521-36 re- 
mains undetected at VHE, a misidentification with the EGRET 
source should be considered. 



6. Discussion and implications on the AGN 
unification scheme 

When applied to AGNs belonging to very diff'erent classes, the 
multi-blob model deduces very similar properties for the size 
of the TeV emitting zones and the values of the magnetic field. 
Furthermore the inferred bulk Lorentz factors F;, can usually re- 
main below a value of 10, thus reconciling SSC models with 
(GR)MHD models, except for the interpretation of the y-ray 
bump of Cen A in terms of synchrotron emission (see the green 
lines in Fig.|7J. A rather unified picture seems therefore to come 
out from the analysis. 

In our scenario, we locate the X- and y-ray emitting region 
in the jet formation zone, with opening angle larger than in the 
global VLBI radio zone, which is more distant from the core 
and mainly located in a region where the jet is more collimated. 
This is somewhat reminiscent of a proposal made several years 
ago by Celotti et al. ( 1993) for unifying X-ray and radio selected 
BL Lac objects. ICelotti et akl d 19931) showed that such a picture 
can be coherent with statistics of the BL Lac population. 
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Since our model accounts for the observation with rather low 
bulk Lorentz factors and with large effective opening angles, we 
have found a way out of the proble m of statistics on the nu mber 
of detected TeV sources invoked bv lHenri & Saugel (I2006h . and 
we can reconcile beamed and unbeamed sources. One important 
consequence of our proposal is that bright radio BL Lacs should 
be TeV emitters. This can be tested by further observations. 

Indeed, Padovani (2007) shows that the common histori- 
cal unified LBL/HBL scheme (the so-called "blazar sequence") 
seems to be ruled out by the discovery of "outliers" low-power 
LBL and high-power HBL sources. Hence it does not seem im- 
possible that objects that are very different at first sight, like 
radiogalaxies such as CenA, could be (faint) VHE emitters. 
Furthermore the fact that M 87 has already been detected in the 
TeV range is encouraging for a future detection of such radio- 
galaxies. In the case of 3C273, this object would rather be a 
misaligned LBL-like object. However one should stress that the 
fact that mainly HBL have been detected up to now at VHE is 
certainly only a selection effect. Thus a TeV detection of 3C 273 
would not be very surprising. 

The scenario we propose here has an additional interest in 
the sense that it allows to solve the long standing paradox of 
the apparent absence of high superluminal motion at the base 
of radio jets of TeV BLLacs. In our model, some X- and y-ray 
emitting plasma blobs are moving close to the line of sight, thus 
allowing for instance Cen A and PKS 052 1 -36 to be potentially 
seen at VHE, while standing in misaligned extended jets. 

This statement a g rees well wi th the most recent studies by 
iGopal- Krishna et al] (l2006l l2007l) who show that viewing an- 
gles, opening angles and Lorentz factors of (sub-)parsec scale 
jets evaluated by radio observations are usually underestimated, 
thus reconciling absent superluminal motions inferred from ra- 
dio observations with high Lorent z facto rs required by different 
families of TeV emission models. IWiital (^006) also underlines 
the possibility to reconcile various contradictory observations by 
considerin g jets with opening angle s of a few degrees. Early 
studies by iDermer & Gehrelsl ([T995h on y-ray observations of 
AGNs with EGRET also indicate that strong beaming is not re- 
quired to account for TeV blazars observations. 

Two different variability time scales appear within our mutli- 
blob scenario. The short time scale is related to the characteris- 
tic size of individual blobs, as already discussed. A longer time 
scale occurs in the case of rotating jets with helicoidal magnetic 
field. This induces a rotation of the cap, and the lag between 
"on-blob" and "inter-blob" emission corresponds to a long term 
y-ray variability. From MHD models, we infer a characteristic 
rotation time of the magnetic field of the order of one year in 
the observer frame for M 87. This modulation could explain the 
variation of the VHE emission observed between low and high 
states in 2004 and 2005. This would imply some periodicity in 
the TeV emission of AGNs, but on timescales hitherto not ex- 
plored. 

7. Conclusion 

We have presented a SSC model to interpret VHE emission 
of M87 as well as other misaligned sources with extended 
optical/X-ray jet. This model accounts in a simple way for a dif- 
ferential Doppler boosting effect by modeling the emission of 
several blobs of plasma located in the broadened formation zone 
of the jet close to the SMBH, just beyond the Alfven surface 
predicted by MHD models. 

Our scenario provides a reasonable interpretation of the 
H.E.S.S. VHE observations of M87 and provides the possibil- 



ity to extend standard leptonic models of TeV blazars to other 
types of AGNs. However, we do not exclude other leptonic or 
hadronic models. For instance, iNeron ov & Aharonian (2007) re- 
cently interpreted the H.E.S.S. observations from 2005 of M87 
by invoking acceleration and radiation of electrons in the black 
hole magnetosphere, which is another kind of leptonic model. 
Hadronic models also cannot be excluded as efficient particle ac- 
celeration processes can occur in the close surroundings of the 
black hole. 

More observations are needed to constrain the emission 
models and especially to distinguish between hadronic and lep- 
tonic scenarii. The upcoming GLA5r mission and the H.E.S.S. II 
project will certainly help to understand the mechanisms at work 
in the AGNs by exploring spectral ranges below TeV, which is 
decisive to constrain the shape of the inverse Compton bump. 
Moreover, the lack of genuine simultaneous multiwavelength 
campaigns on M87 needs to be filled, especially since this 
source is known to be variable at small timescales in VHE. 

Several types of active nuclei are potential emitters of 
VHE photons with predicted TeV fluxes detectable by present 
Cerenkov arrays like H.E.S.S. and MAGIC, or by the next gen- 
eration of instruments such as the CTA project. Such data will 
be crucial to test AGN unifying schemes. 
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